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ABSTRACT 

We have timed four miUisecond pulses, PSRs J1721-2457, J1745-0952, J18I0-2005, and J1918-0642, for up to a total of 10.5 
years each using multiple telescopes in the European Pulsar Timing Array network: the Westerbork Synthesis Radio Telescope in 
The Netherlands, the Nanjay Radio Telescope in France and the Lovell telescope at Jodrell Bank in the UK. The long time span has 
enabled us to measure the proper motions of J1745-0952 and J19I8-0642, indicating that they have transverse velocities of 200(50) 
and 54(7)kms"'respectively. We have obtained upper limits on the proper motion of J1721-2457 and J18I0-2005, which imply 
that they have transverse velocities less than 140 and 400 km s"' respectively. In all cases, the velocities lie in the range typical of 
millisecond pulsars. We present pulse profiles for each pulsar taken from observations at multiple frequencies in the range of 350 to 
2600 MHz, and show that J1810-2005 shows significant profile evolution in this range. Using our multi-frequency observations, we 
measured the spectral indices for all four pulsars, and for J1810-2005 it appears to be very flat. The flux density of J19 18-0642 shows 
extensive modulation which we attribute to the combined effects of refractive and diffractive scintillation. We discuss the possible use 
of including J1721-2457 or J1918-0642 in a pulsar timing array, and find that J1918-0642 will be useful to include when the timing 
precision of this pulsar is improved over the next few years. We have searched archival optical observations to detect companions of 
the binary pulsars, but none were detected. However, we provide lower limits on the masses of the white dwarf companions of PSRs 
J1745-0952 and J1918-0642. 

Key words, stars: neutron - pulsars: general - pulsars: individual: J1721-2457, J1745-0952, J1810-2005, J1918-0642 



1. Introduction 

In this paper we present improved timing solutions for four 
millisecond pulsars (MSPs): PS Rs J 172 1-245 7, J1745-0952, 
and J1918-0642, discovered by 



Edwards & Bailes (2001 ), and 



J18 10-2005, that was discovered byl Camilo etal.| ( |2001 j ). Long- 
term timing of MSPs is an important tool to determine masses of 
the individual stars in binary systems, and by constraining secu- 
lar effects like proper motion, it can be used to improve statistics 
on the evolution of these systems. 

In general, masses of pulsars are not easy to determine. In 
some cases, when the pulsar is in a binary with another compact 
object, high precision pulsar timing observations on extended 
timescales can allow for detecting post-Keplerian parameters of 
the systems, which can be used to separately measure the indi- 
vidual masses of the stars. Space velocities derived from proper 
motion measurements of radio pulsars give clues about the evo- 
lution of these systems and their birth supernovae (e.g.|ToscanOj 
|etaLl ( fT999l l; |Hobbs et al.| ( |2005] l; |Lorimer| ( |2008l l). It is believed 
that recycled pulsars have lower space velocities than normal, 
slowly rotating pulsars. However as MSPs are generally the most 
stable rotators, the effects of their space velocities on their rota- 
tional and orbital parameters will be easier to derive from timing 



measurements of those systems and could be used to deduce the 
intrinsic properties of the stars. 

Another interesting and important use of timing MSPs to 
high precision, is the formation of a pulsar timing array (PTA, 
e.g. Hobbs et al. 2008; Janssen et al. 2008b). An instrument like 
this will use an array of MSPs as the endpoints of a Galaxy-scale 
gravitational wave (GW) detector Current estimates predict that 
to detect a GW background, long-term high precision timing of 



about 20 MSPs is required ( Jenet et al. 2005 ; Van Haasteren et al. 
2008). Increasing the number of stable MSPs in the array will 
improve the detection significance. In order to better understand 
these systems in general, and to determine their suitability for in- 
clusion in a timing array, their long-term timing behaviour needs 
to be determined. Moreover, the best frequency at which they 
should be timed needs to be ascertained. This is a combination 
of the pulse shape at these frequencies and their intensity. We ob- 
served all four pulsars at additional frequecies to find their best 
possible observing frequency for timing purposes, and discuss 
their suitability for inclusion in a PTA. 

Three of the pulsars that are presented in this paper are 
in low-eccentricity binary systems with white dwarf compan- 
ions (see Table[T]l. These systems are usually classified as low- 
mass binary pulsars (LMBPs). The LMBPs distinguish them- 
selves from intermediate-mass binary pulsars (IMBPs) in hav- 
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ing shorter periods (< 10 ms), very low eccentricities (< 10"^), 
and they follow relationships between their orbital period and 
the eccentricity of their orbit, and their orbital period and the 
companion mass ( |Thorsett & Chakrabarty|1999| l. The origin of 
the difference between these two classes of binary MSPs is be- 
lieved to be an evolutionary effect, and mainly due to the pro- 
genitor masses of the companion stars. There are now ~65 MSP 
binary pulsars with probable white dwarf companions known 
(e.g. Lorimer 200 8]) and about 16 of these can be regarded as 
IMBP candidates ( [Jacoby et al.|2006| l 



Detecting an optical counterpart to one of the binary pul- 
sars allows to derive properties of their white dwarf companions. 
Because of differences in the cooling properties of white dwarfs 
in LMBPs and IMBPs, optical observations can in some cases be 
used to distinguish between the two classes (e.g. van Kerkwijk & 
[Kulkarni 1995 ). So far, for all three binary pulsars described in 
this paper, the most recently published limit on optical magni- 
tude of the companion is R > 24 ( van Kerkwijk et al.|2005^ . 



2. Pulsar timing observations and data analysis 

2.1. Westerbork 

The pulsars were observed approximately monthly using the 
pulsar machines at the Westerbork Synthesis Radio Telescope 
(WSRT): PuMa (|Voute et al.|[2002l ) since August 1999, and 
PuMall ( [Karuppusamy et al.|[2008[ ) since August 2006. The 
WSRT has three frequency ranges allocated for pulsar tim- 
ing; for observations centred at 350 MHz we use two bands of 
10 MHz wide, and observations centred at 1380 or 2300 MHz 
use 80 MHz of bandwidth, spread in 8 steps of 10 MHz over 
a range of 160 MHz. The new pulsar machine uses the full 
80 MHz bandwidth at the low frequencies, and all of the avail- 
able 160 MHz bandwidth at 1380 and 2300 MHz. For all four 
pulsars presented here, we mainly used the 1380 MHz band as 
the pulsars were detected in this band, and there are no other re- 
ports of detections in other bands. Also this is the best available 
band for timing at the WSRT as it has the best sensitivity. For 
completeness, we observed all pulsars for one hour at 350 and 
2300 MHz, seeFig[T] 

For a complete description of the PuMa data analysis, we re- 
fer to Janssen et al. (2008a'). Analysis of the PuMall data was 
carried out using the PSRCHIVE ( |Hotan et al.,200 4) software 
suite. For both instruments, each profile was cross-correlated 
with a standard profile at the corresponding observing fre- 
quency (Fig.[T]i, obtained from the summation of high signal- 
to-noise (S/N) profiles, to calculate a time of arrival (TOA) for 
each observation. These were referred to local time using time 
stamps from a H-maser at WSRT and converted to UTC us- 
ing global positioning system (GPS) maser offset values mea- 
sured at the observatory, and GPS to UTC corrections from the 
Bureau International des Poids et Mesures (BIPMQ The TOAs 
were converted to the Solar system barycentre using the JPL 
eph emeris DE40^ We used the tempo2 timing software pack- 
age ( Hobbs et al.|2006 ) to analyse the data. 

Flux densities were calculated for the WSRT observations 
by using the aforementioned bandwidth and observing times, 
the measured S/N ratio for each profile, the gain of the WSRT 
(1.2K/Jy), the pulse duty cycle (~10%), the system temperature 
and the pulsar specfic radiometer equation ( Dewey et al.|1985 i. 
The total system temperature was determined from the synthesis 
data that is taken in parallel with the pulsar observations. 

' http://www.bipm.org 

^ ftp://ssd.jpl.nasa.gov/pub/eph/export/DE405/de405iom.ps 



2.2. Nangay 

PSR J 172 1-2457, J 1745-0952 and J191 8-0642 were observed 
roughly every 3 to 4 weeks with the Nangay Radio Telescope 
(NRT) since 2006. 

Equivalent to a 93-m dish, the Nangay Radio Telescope and 
the BON (Berkeley-Orleans-Nangay) coherent dedispersor were 
used for typical integration times of 45 min. Coherent dedis- 
persion of a 64 MHz band centred on 1398 MHz was carried 
out on sixteen 4 MHz channels using a PC-cluster in the pe- 
riod covered by the observations. The Nangay data are written 
with timestamps which are directly tied to UTC(GPS) in real- 
time, using a Thunderbolt receiver (Trimble Inc.), when convert- 
ing from an analogue to digital signal is performed. Differences 
between UTC and UTC(GPS) are less than 10 ns at maximum, 
and therefore no laboratory clock corrections are needed. One 
TOA was calculated from a cross-correlation with a pulse tem- 
plate for each observation of ~ 45 min. 

The flux densities for the Nangay data were calculated in a 
similar way to those for the WSRT except that the system tem- 
perature was determined from the known receiver temperature 
and sky temperature in the direction of the source. The gain of 
the NRT is 1.4K/Jy. 

2.3. Jodrell Bank 

PSR J18 10-2005 was observed at Jodrell Bank since November 
1997, and PSR J 19 18-0642 since February 2008. Both pulsars 
were observed every two weeks at a centre frequency around 
1400 MHz. The gain of the Lovell telescope is 1 K/Jy. 

The polarized signals from the receiver were were fed into 
an analogue filterbank system with 2 x 32 x 1 MHz channels, 
and filtered and digitised at appropriate sampling intervals and 
incoherently dedispersed in hardware. The resulting dedispersed 
timeseries were folded on-line with the topocentric pulsar period 
and finally written to disc. In the off-line reduction, the two po- 
larisations were summed to form total-intensity profiles. A stan- 
dard high S/N pulse template was used to determine a TOA for 
each observation. During this process, TOAs were referred to the 
local H-maser time-standard and akeady corrected to UTC using 
information obtained via the GPS . 



3. Pulsar timing results and discussion 

We present new timing solutions for the pulsars in Table[T] based 
on 7.5 to 10.5 years of radio pulsar timing as described in Sect. [2] 
The vast majority of our timing data were taken using a cen- 
tre frequency around 1380 MHz, and therefore originally only 
the timing variation across the observing bandwidth was used 
to measure the dispersion measure (DM) of the four pulsars. As 
the combined data sets of WSRT and Nangay or Jodrell Bank 
together cover more than 200 MHz of bandwidth, the DM was 
already determined with good accuracy using only that data. 
However, to allow for refining the DM value even more, and to 
study the pulse shape and flux densities at other frequencies, we 
also observed the pulsars at 350 MHz and 2.3 GHz. This resulted 
in a few detections, see Fig. [T] 

For PSRs J1745-0952, J1810-2005 and J1918-0642, all in 
very low-eccentricity binaries, we have used the ELLl timing 
model (Lange et al. 2001 ) to avoid the high degree of covariance 
between the epoch and longitude of periastron. 

We have measured significant proper motions for 
PSRsJ1918-0642: /^t = 9.2(2) mas yr"' and J1745-0952: 
/jt = 24(2) mas yr"'. By combining those proper motions with 
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Table 1. Timing solutions for the pulsars. 



Pulsar name 


J1721-2457 


J1745-0952 


Jl 8 10-2005 


J19 18-0642 


Fit and data-set 


MJD range 

Number of TO As 

Rms timing residual (pis) 

Weighted fit 

Reduced value" 


52076-54799 
120 
31.9 
Y 

1.05 


52076—54776 
113 
20.8 
Y 

1.02 


50757-54796 
491 
270.7 
Y 

1.01 


52094-54814 
152 
2.6 
Y 

1.04 



Measured Quantities 



Right ascension, a (J2000) 


17h21'"05':4963(3) 


17H5'"09; 1348(2) 


18''10™58^9919(6) 


19higm4gso35270(l3) 


Declination, S (J2000) 


-24°57'06'.'36(6) 


-09°52'39'.'682(15) 


-20°05'08'.'27(14) 


-06°42'34'.'8636(6) 


Pulse frequency, v (s"') 


285.989343507338(15) 


51.609431233273(5) 


30.467142155106(7) 


130.7895141841235(6) 


First derivative of pulse frequency, v (s^^) 


-4.533(5)xl0-'' 


-2.4627(13)xl0-"^ 


-1.3684(ll)xl0-'* 


-4.39525(17)xl0-'* 


Dispersion measure, DM (cm^^pc) 


48.34(3) 


64.27(9) 


241.0(3) 


26.554(10) 


Proper motion in RA, /la (mas yr"' ) 


1.8(1.8) 


-21.2(1.1) 


0(2) 


-7.20(10) 


Proper motion in DEC, fis (mas yr"') .... 


-14(25) 


11(5) 


17(37) 


-5.7(3) 


Binary 




ELLl 


ELLl 


ELLl 


Orbital period, (d) 




4.943453386(12) 


15.01201911(4) 


10.9131777486(12) 


Projected semi-major axis of orbit, x (It-s) 




2.378599(5) 


11.977880(18) 


8.3504716(4) 


TASC (MJD) 




53198.621445(3) 


53195.528458(4) 


53402.53135510(8) 


EPSl (ei) 




0.000009(3) 


0.000009(3) 


-0.00001244(9) 


EPS2 (e,) 




-0.000004(4) 


-0.000017(3) 


-0.00001555(9) 


Set Quantities 


Epoch of frequency determination (MJD) 


53400.0 


53200.0 


53200.0 


53400.0 


Epoch of position determination (MJD) . . 


53400.0 


53200.0 


53200.0 


53400.0 


Epoch of DM determination (MJD) 


53400.0 


53200.0 


53200.0 


53400.0 


Derived Quantities 


Orbital eccentricity, e = yje-^ + 




1.0(0.4)xl0-' 


1.9(4)xl0-^ 


1.991(13)xl0-' 


Omega, a) - arctan (^i/ei) (°) 




114(28) 


152(12) 


218.6(4) 


Characteristic age, (Gyr) 


10.0 


3.3 


3.5 


4.7 


Surface magnetic field strength, (10'*G) . 


1.4 


13.5 


22.4 


4.4 


Distance*, (kpc) 


1.3(2) 


1.8(3) 


4.0(4) 


1.23(13) 


Spin period (ms) 


3.49663 


19.3763 


32.8222 


7.64587 


Spin period derivative 


5.55x10-21 


9.23x10-20 


1.47x10-" 


2.57x10-20 


Mass function, (Mq ) 




0.000591270(4) 


0.00818737(4) 


0.0052494438(8) 


Minimum companion mass, (Mq ) 




0.11 


0.28 


0.24 


Total proper motion, fij (mas yr"' ) 


15(25) 


24(2) 


17(37) 


9.2(2) 


Transverse velocity, Vx (kms"') 


< 140'' 


200(50) 


< 400^ 


54(7) 



Assumptions 


Clock correction procedure 


TT(TAI) 


Solar system ephemeris model 


DE405 


Model version number 


5.00 



Note: Figures in parentheses are the nominal Icr tempo2 uncertainties in the least-significant digits quoted. 

"Before combining the data sets from the dilferent observatories, the errors on the TOAs were scaled with a constant factor to have each individual 
data set return a reduced ~ I - ''The DM distances are estimated from the Cordes & Lazio ( 2002 i model. ^ Velocity limits from Eq.fT] 



the DM-derived distances from the |Cordes & Lazio|2002| model, 
we have calculated transverse velocities for those systems. For 
the other two pulsars, we use the marginally significant values 
from the timing solution to calculate a limit on their transverse 
velocities, see Sect.|3.1| 



3.1. Proper motion and velocity 

We use our proper motion measurements to derive transverse ve- 
locities, or limits on those, for all four pulsars (see Table[T]i. We 
have a highly significant measurement of the proper motion of 
PSR J 191 8-0642, and a significant detection in right ascension 
for PSR J 1745 -0952. The proper motion values as presented 
for the other two pulsars are not significant. We note that both 
PSR J 172 1-2457 and J 18 10-2005 have ecliptic latitudes very 
close to the ecliptic plane (J3 = -1°8 and 3^3 respectively), and 



therefore their proper motions may be better and more indepen- 
dently fit in ecliptic parameters. However, as the coordinate sys- 
tems are nearly parallel because both pulsars are near a - 18^, 
fits in this coordinate system did not result in an improvement 
of the limits on their proper motions, or the significance of their 
positions. As both of these pulsars lie quite close to the ecliptic, 
it may be that they will show effects due to the solar wind (You 
et al.|2007 1. However at the present timing precision we see no 



influence on the observed TOAs. 

Due to changes in the projected line of sight, transverse 
motions of pulsars can affect observed peiiodicities in the sys- 
tem ( Shklovskii| 1 97Q\ . We present the contributions to the spin 
period derivative of each pulsar in Table|2] For completeness, 
the contributions to the observed P from accelerations in the 
Galactic potential are also presented, although those are not im- 
portant effects for these pulsars. To calculate the magnitude of 
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Table 2. Contributions to the period derivative. Contributions 
due to accelerations in the Galactic potential, Pcnh are calculated 
using equations in|Toscano et al. ( 1999 1. 



PSR 




-Pobs 




^Shk.max 


P Gal,± 




P Gain 




J1721 


-2457 


5.55x10- 


-21 


2.5 X 10- 


21 


4.0x10- 


23 


4.4x10- 


-22 


J 1745 


-0952 


9.23x10- 


-20 


4.9 X 10- 


20 


4.5x10- 


22 


2.8x10- 


-21 


J1810 


-2005 


1.47x10- 


-19 


5.4 X 10- 


20 


1.1x10- 


23 


1.8x10- 


-20 


J1918 


-0642 


2.57x10- 


-20 


1.9 X 10- 


21 


1.4x10- 


22 


3.7x10- 


-22 



the Shklovskii term we have used the total proper motion values 
as presented in Table[T] For PSR J1918-0642, the contribution 
to its spin period derivatives, Pshk, from proper motion is less 
than 8% of the observed value. However for PSR J 1745 -0952, 
the effect can be as large as 50% of the observed value, leading 
to an underestimation of the characteristic age and an overesti- 
mation of the magnetic field. Moreover, as the proper motions 
for the solitary pulsar J 172 1-2457 and PSR J1810-2005 have 
not been determined significantly yet, the Shlovskii term may 
even completely dominate the observed P for those pulsars. To 
calculate an upper limit for the proper motion, and accordingly 
a limit to the transverse velocity, we can use the observed P as a 
maximum: 



(1) 



This yields fij < 23masyr ' and Vj < 140 km s 'for 
PSR J 172 1-2457 and jUt < 22masyr-i and Vt < 400kms-i 
for PSR J 18 10-2005. 

About 50% of the known solitary and low-eccentricity bi- 
nary pulsars have measured transverse velocities (Tables 2 and 
4 in [Lor imer 2008). [Hobbs et aL ( 2005[ l quote mean values of 
2D speeds of 77(16) kms^' for solitaiy MSPs and 89(15) kms^' 
for binary MSPs. Compared to those, the (limits on) transverse 
velocities that we have derived from the proper motion measure- 
ments represent normal velocities for recycled pulsars. 













< 








Shk 




Obs 



3.2. Profiles 

Apart from our normal timing observing frequency of 
1380 MHz, we have observed all four pulsars at additional fre- 
quencies: 350 MHz and 2300 MHz. The results are shown in 
Fig. [T] It is now standard procedure for WSRT timing obser- 
vations of MSPs to use both pulsar machines, and as we have 
twice the bandwidth available in PuMall, we present the pro- 
file of the detected pulsars at 350 MHz and 2.3 GHz from that 
data, see Fig. [T] Single observations of 55 min. were used to 
generate the profiles at 350 and 2300 MHz. Where no profile is 
plotted, the pulsar was not seen at that frequency. The profiles at 
1380 MHz were generated using one year of regular timing ob- 
servations and consist of 10 to 12 observations of 25 min, except 
for PSR J1918-0642, where again one observation of 55 min. 
was used. PSR J 1745-0952 was detected with NRT at 2600 GHz 
in an observation of 50 min. 

As reported by|Kramer et al. (1998]), for most MSPs there is 
very little development of pulse profiles with frequency, which 
can be explained by MSPs having a very compact magneto- 
sphere. For all pulsars except PSR J 18 10-2005, we indeed see 
no changes in the pulse profile with changing frequency, ex- 
cept for changes in the ratios of the heights of some compo- 
nents, within the significance of our detections. Compared to 
its profile at 1380 MHz, the 2300 MHz profile of 1810 ap- 



pears to have more components in the main peak. Also the pro- 
file appears to show a leading component about 70 degrees be- 
fore the main peak, although higher S/N will be needed to con- 
firm this. Of all the MSP profiles presented by [Kramer et al.| 
(1999 1 none show more complex profiles at higher frequencies 



than at lower frequencies. Thus a further study of the profile 
of PSR J1810-2005 and its polarisation properties could reveal 
more about MSP emission profiles. The extra features in the peak 
of the high-frequency profiles of PSR J18 10-2005 could be an 
indication that the profile at 1380 MHz suffers from scattering 
effects. However, the Cordes & Lazio (2002) model predicts very 



little scattering for this pulsar at both 1380 MHz (0.26 ms) and 
2300 MHz (0.03 ms) suggesting that the features in the main 
peak of the pulse profile are intrinsic to the pulsar It may turn out 
that, when better sensitivity can be achieved at higher frequen- 
cies, PSR 1810-2005 will give better timing results compared 
with the more commonly used 1380 MHz as its profile shows 
more features at higher frequencies. 

3.2.1 . Useful for pulsar timing array? 

PSR J 172 1-2457 shows very stable rotational behaviour All pa- 
rameters of the timing solution presented in Table T] are consis- 
tent with the original timing solution as published by [Edwards 
j& Bailes] ( [200 l| l. However, for the present S/N ratio, the wide 
profile of this pulsar does not allow for the high timing precision 
that is required for PTA pulsars ( Jenet et al.|2005) l. 

In contrast, despite its relatively long spin period compared 
to most pulsars considered for PTA studies, PSR J 191 8-0642 
may be worthwile to include. Already at 1380 MHz the profile 
has a very sharp peak, allowing for very precise timing over 7.5 
years (rms - 2.3 yus). Furthermore, it is possible that better re- 
sults for pulsar timing arrays may be deduced by including not 
only the best-timing pulsars, but extending the PTA pulsar set 
with pulsars that have timing solution with rms ~ 1 fis { Jenet 
et al. 2005 Van Haasteren et al. 2008). The pulse profile as 
shown in FigT|T]was generated by the PuMall machine in WSRT, 
promising even better timing precision in the ongoing timing 
programme. Observations with PuMall use twice the bandwidth 
of those with PuMal and so give an increase in the sensitivity. 
Moreover, because PuMall uses coherent dedispersion, it yields 
a sharper profile. The combination of these effects should lead 
to an improvement in the timing precision of a factor between 2 
and 4 for this pulsar over the next five years which may bring it 
to the required level for a PTA pulsar 



3.3. Flux density and scintillation 

We calculated flux densities for the four pulsars as described in 
based on the profiles as shown in Fig.[T| see Table[3] 



Sect. 2.1 



We note that apart from those at 21 cm, the profiles are gener- 
ated from single observations only, and therefore the resulting 
flux densities should not be regarded as accurate measurements, 
but are indicative of the intensities of the pulsars at the observed 
frequencies. The spectral indices for these pulsars are normal ex- 
cept for PSR J1810-2005, which appears to have a relatively flat 
spectrum. We observed this pulsar a couple of times at 2.3 GHz, 
and the resulting flux densities were varying by a factor 2, in- 
dicating that the brightness of this pulsar may be affected by 
scintillation. 

Not long after we began our observations of 
PSR J1918-0642 we noticed that there were a number of 
occasions when we either did not detect the pulsar at all or it 
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Pulse phase 



Pulse phase 



Pulse phase 



Pulse phase 



Fig. 1. Profiles for the PSRs at multiple frequencies. All profiles were generated from a single 55 min. observation using PuMall, 
except the 1380 MHz profiles of PSRs J1721-2457, J1745-0952 andJ1810-2005 that used 10 to 12 hours of PuMall observations, 
and the 2600 MHz profile of PSRJ1745-0952 was generated from a 50 min. Nan§ay observation. These profiles are to show their 
appearance; for timing purposes we used high S/N templates as described in the text. 



was extremely bright. A plot of the flux densities determined 
from all observations is shown in Fig.|2] 

The left hand panel shows a histogram of flux values from 
both Nanfay and WSRT observations with the non-detections in 
the WSRT data shown in the far left hand bin. The right hand 
panel shows the variation in flux density as a function of time 
indicating the variations are typically on short timescales, al- 
though a long period of non-detections can be seen in 2006- 
2007. We note that the pulsar is not detected in approximately 
20% of all WSRT observations. The pulsar is always detected, 
although sometimes with very low S/N ratio, in the Nangay data. 
This is not suprising due to the higher gain and longer integra- 
tion times for the Nan9ay observations meaning that the NRT is 
approximately 30% more sensitive than the WSRT for this pul- 
sar. 

If these flux variations are due to scintillation then the 
quite frequent non-detections would suggest that the scintilla- 
tion bandwidth must be at least as wide as the PuMal bandwidth 
of 80 MHz. For the dispersion measure of this pulsar at this ob- 
serving frequency such a large scintillation bandwidth seemed 

t|2001| l. Checking the literature. 



too large to be diffractive ( Rickett 



we found four pulsars with similar dispersion measures which 
had measured diffractive scintillation bandwidths in the range 
15 - 70 MHz ( [Johnston et al . 1998 , Wang et al.|2005[ ) at 21 cm. 
While this is a large range the diffractive scintillation bandwidth 



would need to be above the highest measured value to explain 
the observed flux variations. 

Using the improved bandwidth of our PuMall measure- 
ments we were able to better explain the scintillation proper- 
ties of PSR J1918-0642. In Fig.|3]we show the dynamic spec- 
tra from two observations of PSR J 19 18-0642, separated by just 
one month. The left hand side shows a classic case of scintil- 
lation with a scintillation bandwidth of approximately 5 MHz. 
Comparing this value to the range of bandwidths seen for sim- 
ilar dispersion measure pulsars, as discussed earlier, it is a fac- 
tor of three smaller than the lowest measured scintillation band- 
width. In contrast, the right hand plot shows just a single scintle, 
which is perhaps somewhat broader than 5 MHz, across the full 
160 MHz of bandwidth. It may be possible that these variations 
are purely statistical in nature, although the lack of any scin- 
tiles in the band of at least 80 MHz in width, appears to happen 
about 20% of the time. We note that |Gupta et al.| ( |I994"l l explain 
variations in the scintillation properties like this as being to the 
refractive modulation of the diffractive interstellar scintillation 
pattern. While beyond the scope of this paper, the techniques 
1994[ l could be used to further exam- 



outlined by Gupta et al. 



ine the relationship between the velocities in this system and the 
interstellar medium along the line of sight. We note also that 
very wide bandwidth observations of this pulsar provide the best 
opportunity to do high precision timing, as when the pulsar is 
bright it can be timed to high precision. 
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Fig. 2. Flux density variations of PSR J1918-0642 obtained with the WSRT and NRT from observations made at 21 cm. The left 
panel shows the distribution of flux densities from the two observatories, where the leftmost column represent all non-detections. 
In the right panel we show the individual flux density values as a function of time, where the length of the line in each case is 
proportional to the flux density, as indicated in the legend. Observations where there was no clear detection (ND) of the pulsar and 
thus only an upper limit on the flux density can be determined, are plotted with a cross. 




- 20 40 60 - 10 20 30 



Time (minutes) Time (minutes) 

Fig.3. Dynamic spectra of PSR J1918-0642 plotted as a greyscale of intensity as a function of time and frequency. Darker pixels 
correspond to higher intensity values and the intensity scales linearly with the levels of grey. The dynamic spectra are from two 
observations made using PuMall on 23 February 2008 (left) and 24 March 2008 (right). 



Table 3. Flux densities and spectral index. 



PSR 




S350 (mJy) S 


1380 (mJy) S2273 (mJy) 


SI 


J1721- 


-2457 




0.58(2) 0.28(3) 


-1.5(2) 


J1745- 


-0952 


1.8(3) 


0.38(3) 


-1.14(12) 


J1810- 


-2005 




1.33(2) 0.95(4) 


-0.67(3) 


J1918- 


-0642 


5.9(6) 


0.58(2) 


-1.67(6) 



Note: These numbers correspond to the WSRT profiles in Fig. [T| 
Measurements of flux densities at 350 and 2273 MHz, and 
PSRJ1918-0642 at 1380MHz are based on one or two 50 min. ob- 
servations, and therefore should be regarded as indicative of the inten- 
sities of the pulsars at the observed frequencies only. The ranges in flux 
densities, and consequently, SI calculations, are now based on the un- 
certainties in duty cycles. 



4. Optical observations 

We have analysed archival observations of PSR J 1745 -0952 and 
PSR J1810-2005. These observations were obtained with the 
ESO Multi-Mode Instrument (EMMI) at the 3.6 m telescope of 
the European Southern Observatory at La Silla, Chile on March 
15th, 2004. The instrument has two 2Kx4K detectors which 
were read out with 2x2 binning, providing a pixel scale of 
0'.'33pix"'. The observations consisted of two 10 min. /?-band 
exposures of PSR J 1745-0952 and two 5 min. 7?-band exposures 
for PSR J18 10-2005. For both pulsars, these long exposures 
were preceded by short (10 s) exposures, also in the /?-band. The 
seeing during these observations was about 0'.'86 to 0'.'88. All 
images were bias-subtracted and flat-fielded using twilight flats. 
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Fig. 4. 45" X 45" subsections of the deep images of PSR J1745-0952 (left) and PSR J1810-2005 (right). To retain visibility, the 
pulsar positions are denoted by a circles with a 1" radius, even though the actual error ellipses are considerably smaller than this. 



The two deep images of each pulsar were averaged and used for 
further analysis. 

The short 10 s observations were astrometrically calibrated 
against the 2nd version of the USNO CCD Astrograph Catalog 
(UCAC2, [Zacharias et al.|2004| . For the PSR J1745-0952 obser- 
vations, 34 UCAC2 stars overlapped with the chip on which the 
pulsar was located, of which 26 were not saturated and appeared 
stellar and unblended, giving rms residuals of 0'.'056 in right as- 
cension and (J.' 091 in declination. A secondary astro me trie cata- 
log was created from the stars on the short exposure and used to 
calibrate the average of the two 10 min. exposures. For the trans- 
fer, about 300 stars were used, giving residuals of 0'.'022 in right 
ascension and 0'.'025 in declination. Combined with the uncer- 
tainty on the timing position of the pulsar, the quadratic uncer- 
tainty of the pulsar position on the deep optical image is 0'.'060 in 
right ascension and ^.'096 in decUnation. For PSR J1810-2005, 
we used 23 of the 32 overlapping UCAC2 stars, giving rms resid- 
uals of 0'.'049 in right ascension and 0'.'091 in declination. The 
transfer of the secondary catalog to the deeper average of the 
two 5 min. 180 secondary standards, giving residuals of 0'.'028 in 
right ascension and 0'.'018 in declination. The final uncertainty 
on the pulsar position on the deep image is 0'.'058 in right as- 
cension and 0'.'220 in declination. No sources are present in the 
error circles of both pulsars, as shown in Fig.|4] 

In order to determine upper limits on the brightness of the 
pulsar companions, we performed PSF photometry on the av- 



eraged deep images using DAOPHOTII (Stetson 1987) 1. The 
resulting instrumental magnitudes were calibrated against 47 
standard stars in the PG 1525-071 using the calibrated magni- 
tudes of Stetson ( 2000|. Since only /?-band observations were 



taken, only the zeropoint was fitted for the calibration, giving 
rms residuals of 0.02 mag. The standard La Silla 7?-band ex- 
tinction of 0.07 mag per airmass was used to correct for extinc- 
tion between the standard observations taken at an airmass of 
1.08 and the pulsar observations taken at an airmass of 1.34 (for 
PSR J 1 745 -0952) and 1 .25 (for PSR J 1 8 10-2005). Based on the 
magnitudes and uncertainties of faint stars in the images, we es- 



timate the 3cr detection hmit at 7? > 24.49 for PSR J1745-0952 
and R > 24.02 for PSR J1810-2005. 



4. 1 . Model predictions 

Though no optical counterparts to PSR J1745-0952, 
PSRJ1810-2005 and P SRJ1918-0642 (R > 24, see 
van Kerkwijk et aT]|200 5) were detected, it is still worthwhile 
to compare the upper limits on the white dwarf magnitudes 
with those predicted by models. To convert the upper limits on 
the apparent magnitudes to limits on the absolute magnitude, 
we need estimates for the distance and absorption. Estimates 
for the distance to each pulsar using the observed dispersion 
measure and the |Cordes & Lazio| ( |2002[ ) model for the Galactic 
distribution of electrons are given in Table[T] The V-band 
absorption Ay along the line-of-sight and at the distance of each 
pulsar was estimated using the Dri mmel et al.| ( |2003 1 model 
for Galactic extinction and converted to /J-band absorption 
As using the extinction coefficients of Schlegel et al. ( 1998| l 
(Ar = 0.819Av). We estimate upperlimits on the absolute ^- 



band magnitude of Mr > 11.3 for PSR J1745-0952, Mr > 3.3 
for PSR J18 10-2005 and Ms > 12.7 for PSR J1918-0642. 
In Fig.|5] we show predictions of absolute magnitude and 



cooling ages from the white dwarf cooling tracks by Rohrmann 
|et al. ( 2002[ i and Bergeron et al. ( 1995[ l. The helium-core white 



dwarf cooling tracks show a dichotomy in the cooling proper- 
ties of helium-core white dwarfs caused by differences in the 
thickness of the hydrogen layer surrounding the helium core of 
these white dwarfs (see, e.g. |Sch5nberner et al.| 2000[ |Althaus| 
et al. 2001). White dwarfs heavier than 0.18-0.20 Mq have thin 
envelopes and cool faster than lighter white dwarfs which have 
thick envelopes. Cooling tracks appropriate for heavier carbon- 
oxygen white dwarfs with hydrogen-rich atmospheres are even 
fainter 

Under the assumption that the characteristic pulsar ages are 
equal to the white dwarf cooling ages, we plot the upper limits 
on our white dwarf pulsar companions and detections of others 
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Cooling oge (Gyr) 

Fig. 5. Absolute 7?-band magnitude and cooling age predic- 
tions. Shown are helium-core white dwarf cooling tracks by 
Rohrmann et al.| ( 2002 1 and carbon-oxygen white dwarf cool- 
ing tracks with hydrogen-rich atmospheres by Bergeron et al. 
( |1995 |l (log g - 8 corresponds to ~ 0.6 Mq while log g = 9 
corresponds to ~ 1.2 M©). Helium-core models with masses 
below approximately 0.1 8-0.20 Mq have thick hydrogen en- 
velopes and continue to residually burn hydrogen, keeping the 
white dwarf hot and slowing the cooling. Heavier models have 
thin hydrogen envelopes and cool much faster, creating a di- 
chotomy in the cooling properties of white dwarf companions 
to millisecond pulsars. White dwarfs with carbon-oxygen cores 
cool even faster. Shown with error bars are PSRJ02 18-1-4232 
( |Bassa et al.||2003Tl, PSRJ0437-4715 ( |Danziger et aLl[T993) , 
] >SRJ0751-Hl807 passa eFaLllIOOel and PSR J1012-I-5307 
( Lorimer et al.|1995[ l, where the white dwarf cooling age is as- 
sumed to be equal to the characteristic age of the pulsar (us- 
ing braking indices of n - 3). Shown with upper limits on the 
absolute /?-band magnitude are the values for PSR J 1745 -0952 
and PSR J1918-0642. For these two pulsars, thick hydrogen en- 
velopes are excluded, indicating the white dwarf companions in 
these systems are heavier than about 0.2 Mq. 



in Fig.|5] Because of the large distance and absorption towards 
PSR J1810-2005, the limit does not constrain any white dwarf 
parameters. For PSR J1745-0952 and PSR J1918-0642 the up- 
per limits exclude white dwarfs with thick hydrogen envelopes. 
As such, the mass of the white dwarf in these systems is con- 
strained to > 0.2 Mq. For PSR J1745-0952 this mass Hmit, 
combined with the low massfunction of the system, constrains 
the inclination of the binary orbit to the low value of / < 34° 
(assuming Mps, - 1.35 Mq). Should this pulsar have a carbon- 
oxygen white dwarf, which is not excluded by the optical obser- 
vations, the hmit on the inclination will be even lower. 

5. Conclusions 

Using the WSRT, NRT and Lovell telescopes, we have 
timed four MSPs for 7.5 to 10.5 years. We have pre- 



sented updated timing solutions, pulse profiles at multiple 
frequencies for each pulsar, and scintillation parameters for 
PSR J1918-0642. We have measured transverse velocities for 
PSRs J 1745-0952 and J1918-0642, and set Hmits on the veloc- 
ities of PSRs J 172 1-2457 and J18 10-2005. All velocities are 
consistent with previously published distributions for solitary 
and binary MSPs. 

We have analysed archival optical observations for the bi- 
nary MSPs and found no companions to the pulsars. From the 
magnitde limits we deduce for the companions, we can exclude 
white dwarfs with thick atmospheres. This indicates that the 
companions must be heavier than about 0.2 Mq . At this point, 
the mass restrictions as well as the optical magnitude limits 
give no conclusive information to classify PSRs J 1745 -0952 or 
J 18 10-2005 as either LMBPs or IMBPs. 

For low-eccentricity binary pulsars, the only post-Keplerian 
parameter that is likely to be measurable is the Shapiro delay. 
The upper limit for inclination of PSR J 1745 -0952, i<34 , sug- 
gests that the effect of Shapiro delay in timing will be very low 
and therefore we can not expect to use pulsar timing analyis to 
disentangle the individual masses of this system. However, the 
expected luminosity of M« > 11.3 indicates that using a dedi- 
cated optical observation, the companion of this pulsar may be 
detectable, and this will therefore be the most promising way of 
deducing the pulsar and companion masses. 
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